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Abstract By comparing the PSCz galaxy power spectrum with haloes from nested 
and phased N-body simulations, we try to understand how IRAS infrared- 
selected galaxies populate dark-matter haloes. We pay special attention 
to the way we identify haloes in the simulations. 

The galaxy power spectrum measured by Hamilton & Tegmark (2002) 
from the PSCz redshift survey (Saunders et al, 2000) has a remarkable 
power-law form over 4 decades of wavenumber. Interpreting it in the 
ACDM paradigm requires scale-dependent bias between galaxies and 
dark matter, which we hoped to reproduce in N-body simulations. 

We ran four 256 3 -particle AP 3 M dark matter simulations of comov- 
ing boxsize 256, 128, 64, and 32 h~ l Mpc, each with a fixed Eulerian 
softening length of 10 /i -1 kpc. The boxes are nested in the sense that 
the phases of initial fluctuations match in the centers of each box, en- 
suring similar structures there. (For a movie depicting this, look at 
http : //casa. color ado . edu/ ~neyrinck / nesthalf . mpg) 

The non-trivial step in this project is to go from a collection of dark 
matter particles to a set of haloes that could possibly host real galaxies. 
We used DENMAX (Gelb & Bertschinger 1994) with a smoothing length 
of one fifth the mean interparticle separation, which gives a mass spec- 
trum similar to that given by Press-Schechter (1974). We then varied a 
low mass cutoff in the list of haloes to obtain the best fit to PSCz. How- 
ever, there seemed to be an exclusion effect of close pairs because of the 
rather large smoothing length, showing itself as a small-scale downturn 
in the halo power spectra. 

We thus tried "DENMAX 2 ," running DENMAX with half the previ- 
ous smoothing length on each halo separately to resolve subhaloes. We 
wanted to characterize sets of DENMAX and DENMAX 2 haloes with the 
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same parameter (mass would no longer do, since the mass returned by 
DENMAX changes with smoothing length), so we used central density, 
which turns out to be well-correlated with mass in our simulations. 

As expected, DENMAX 2 resolved more close pairs than DENMAX, 
extending the correlation function power law at the small-scale end by 
a factor of two or three. Fig. 1 shows the best fits from each of the 
four nested simulations. The central density cutoff in each simulation 
corresponded to a mass cutoff of about 1O U M . 

In conclusion, we have succeeded in reproducing the scale-dependent 
bias of the PSCz power spectrum at scales ^> 30 h^ 1 kpc with N-body 
simulations. Extending the analysis to the smallest scales requires re- 
solving the substructure of dark matter haloes, which we have done with 
an algorithm we call DENMAX 2 . 
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Figure 1 The PSCz galaxy power spectrum and its errors, along with best fits from 
the 256, 128, 64, and 32 h~ x Mpc simulations (error bars suppressed). Each fit is 
good over a range of ~256 in wavenumber. 



